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The Hamiltonian structure of spacetimes with two commuting Killing vector 
fields is analyzed for the purpose of addressing the various problems of time that 
arise in canonical gravity. Two specific models are considered: (i) cylindrically 
symmetric spacetimes, and (ii) toroidally symmetric spacetimes, which respectively 
involve open and closed universe boundary conditions. For each model canonical 
variables which can be used to identify points of space and instants of time, i.e., 
internally defined spacetime coordinates, are identified. To do this it is necessary 
to extend the usual ADM phase space by a finite number of degrees of freedom. 
Canonical transformations are exhibited that identify each of these models with 
harmonic maps in the parametrized field theory formalism. The identifications 
made between the gravitational models and harmonic map field theories are com- 
pletely gauge invariant, that is, no coordinate conditions are needed. The degree 
to which the problems of time are resolved in these models is discussed. 
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1. Introduction. 



Throughout a large part of the history of general relativity much effort has been ex- 
pended toward disentangling the true degrees of freedom of the gravitational field from 
the "pure gauge" degrees of freedom brought into the theory via the principle of general 
covariance. In the Hamiltonian form of the theory this problem involves understanding 
the solution space of the initial value constraints and the appropriate free data for the 
Cauchy problem. Classically, a characterization of the true degrees of freedom is relevant 
for analyzing dynamical evolution of strongly gravitating systems, e.g., binary black hole 
systems, as well as for understanding fundamental issues in relativity, e.g., cosmic censor- 
ship. The intertwining of gauge degrees of freedom and dynamical degrees of freedom is 
especially vexing in quantum gravity, where it leads to many of the "problems of time" 
[1]. Evidently, canonical quantization of the gravitational field would be expediated by 
a sufficiently explicit characterization of the true degrees of freedom. For the most part, 
the strategy for doing this stems from the original Hamiltonian description of gravitation 
provided by Arnowitt, Deser, Misner [2] and Dirac [3]. The philosophy adopted there is 
that the Einstein field equations define an "already parametrized field theory" in which 
certain non-dynamical canonical variables represent points of space and instants of time, 
relative to which the true degrees of freedom evolve. In the ADM approach, the dynamical 
content of relativity is exposed by coordinate conditions which fix the non-dynamical gauge 
variables. This leads to a desciption of gravitational dynamics relative to a fixed foliation 
of spacetime, that is, relative to a fixed family of observers. A prominent example of this 
approach appears in the conformal approach to the initial value problem [4]. From such 
work it seems that the ADM approach is quite adequate for addressing many problems in 
classical relativity. However, the price paid for obtaining technical control over the true 
degrees of freedom is that general covariance is lost in the sense that one is obliged to 
view dynamics from the point of view of a given set of observers. In quantum theory this 
provides a first instance of a problem of time, which might be called the problem of general 
covariance: How to give the state of the gravitational field on an arbitrary hypersurface, 
that is, with respect to arbitrary observers? In the ADM approach one is prohibited from 
even asking this question. 

An alternative approach to describing the true degrees of freedom of the gravitational 
field that preserves general covariance is available. This approach still relies upon the pos- 
sibility of extracting "many-fingered time" degrees of freedom — or more precisely, spacelike 
embeddings of Cauchy surfaces from the gravitational phase space, but does not fix the 
foliation with coordinate conditions. Instead, one describes evolution of the true degrees of 
freedom relative to an arbitrary foliation, i.e., one casts the Einstein equations in the form 
of a parametrized field theory [5] . This point of view was developed in considerable detail 
by Kuchaf [6] , who called it the "bubble time" dynamics of the gravitational field. Many 
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of the problems of time are mitigated using this approach [1], which is now known as the 
"internal time formalism" . Implementation of the internal time formalism hinges upon the 
possibility of (i) finding a canonical transformation on the gravitational phase space which 
separates four canonical variables to play the role of spacelike embeddings of a Cauchy 
surface into spacetime; (ii) eliminating the momenta conjugate to the embeddings by the 
initial value constraints. If this can be done, the remaining variables represent the true 
degrees of freedom, whose dynamical evolution occurs via correlation with the arbitrary 
spacelike slices provided by the embedding variables. So far, the internal time approach 
to gravitational dynamics has been implemented in lower-dimensional models, typically 
symmetry reductions of the full theory. There are a plethora of homogeneous cosmological 
models in which it is possible, at least locally, to isolate a canonical variable representing 
time and then solve the constratints for the momentum conjugate to time. These models 
possess a finite number of degrees of freedom. Relatively few field theoretic models, i.e., 
models possessing an infinity of degrees of freedom, exist in which one can implement the 
internal time program. Such models are, to our knowledge, always 2-dimensional field 
theories, most notably Einstein- Rosen waves [7] , plane-gravitational waves [8] , spherically 
symmetric vacuum gravity [9], the bosonic string [10,11], and related models [12]. For the 
most part, these models are generally covariant theories of one or more free fields, e.g., the 
dynamics of the Einstein-Rosen wave is that of a single free scalar field representing the 
metric amplitude of the wave. 

Our goal in this paper is to establish an internal time formalism for the reduction of the 
vacuum Einstein equations obtained by assuming the existence of two commuting spacelike 
Killing vector fields. We consider models involving both open universe and closed universe 
boundary conditions. Our work generalizes previous models, in particular the Einstein- 
Rosen waves [7], in two key ways. First, the reduced system of equations describes a pair 
of interacting fields. To our knowledge, this is the only symmetry reduction of the vacuum 
Einstein equations known that admits an internal time formulation and in which the true 
degrees of freedom constitute a non-linear field theory. Second, in terms of the two Killing 
vector model, we illustrate in detail the fact that, strictly speaking, general relativity is not 
an already parametrized field theory. In order to extract canonical variables representing 
embeddings from the ADM phase space it is necessary to extend that phase space by a 
finite number of degrees of freedom. For open universes, the additional non-gravitational 
variables represent the asymptotic location of a spacelike hypersurface relative to an inertial 
frame "at infinity" [13]. For closed universes, the intrinsic and extrinsic geometry of a 
hypersurface are inadequate to define the embedding of that hypersurface in spacetime [14] ; 
additional non-gravitational variables are needed to accomplish this. Given the extensions 
of the phase space, the resulting dynamical system is in each case a generally covariant 
(parametrized) formulation of harmonic maps from a flat 3-dimensional spacetime to a 
2-dimensional target space of constant negative curvature. 
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In the next section we define the models. Our open universe model is obtained by 
imposing cylindrical symmetry on spacetime. Our closed universe model, first considered 
by Gowdy [15], defines the universe to be a 3-torus which admits a 2-torus group of 
isometries. We do not impose hypersurface orthogonality (reflection symmetry) on the 
Killing vectors and the resulting dynamical theory is in each case intrinsically non-linear, 
in contrast to the Einstein- Rosen waves [7] or "polarized" Gowdy models [16]. In the 
language of the linearized theory, both gravitational polarizations propagate and interact. 
In §3 we exhibit canonical transformations, inspired by [7] and [10], which allow us to 
extract the embeddings, solve the constraints, and reveal the true degrees of freedom in 
each model. Other two Killing vector models, which differ in the choice of isometry group, 
spacetime topology, etc. , can be treated in much the same way as we do here. As mentioned 
above, these results lead to substantial simplifications in the problems of time that occur 
in canonical quantum gravity. We discuss these problems and the general structure of the 
quantum theory based on the internal time formulation of the two Killing vector models 
in §4. In the cylindrically symmetric model it is necessary to keep track of the behavior 
of fields on the axis of symmetry and at infinity. We summarize our boundary and fall-off 
conditions in an appendix, which is referred to throughout the paper. 

2. Spacetimes with two commuting Killing vectors. 

We are going to study the Hamiltonian structure of spacetimes (Ai,g) which admit 
two linearly independent spacelike Killing vector fields K^, A = 1,2. The Killing vector 
fields are assumed to commute: 

[K A ,K B ] a = 0, (2.1) 

and to generate a two-dimensional "orthogonally transitive" group Q of isometries. This 
latter requirement means that the distribution of two dimensional vector spaces spanned by 
K ^ at each point is orthogonal to a foliation of M. by surfaces M M. Our assumptions 
amount to demanding that the integral curves of the Killing vector fields provide a fibration 
7r: M — > M of the spacetime manifold by spacelike surfaces. The two-dimensional manifold 
M is the "space of orbits" of the Killing vector fields. 

The spacetime metric g a p defines three functions 

AAB = QapKlK^ (2.2) 

representing the lengths and inner products of the Killing vector fields. Because the Killing 
vectors are spacelike and linearly independent at each point, this symmetric matrix will 
have a symmetric inverse \ AB . We can use \ AB to define a projection operator 7^ into 
the (co-)tangent space to each point of M: 

7? = $ - * AB K%K Bf3 . (2.3) 
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The associated tensor field 



7a/3 = 9a(3 ~ X AB K Aa K Bf} , (2.4) 



satisfies 

£K A -Y a p = and K% laf3 = 0. (2.5) 

There is a one-to-one correspondence between tensor fields on M and tensor fields on M 
with vanishing Lie derivatives along and which are completely "orthogonal" to 
[17,18]. This correspondence means that r ) a p uniquely defines a Lorentz metric on M. 

Coordinates x a = (x A , x a ), A = 1, 2, a = 3, 4 can be chosen such that 

K a A = 5% (2.6) 
and the spacetime metric takes the form 

g a pdx a dx /3 = \ AB dx A dx B + <y ab dx a dx b , (2.7) 

where \ AB = X AB (x a ), and 7 a5 = ~fab(x a )- 

There are only six independent Einstein equations because the orthogonal transitivity 
requirement forces 

^K^ A G af3 = 0. (2.8) 
Thus the non-trivial Einstein equations can be taken to be 

K a A K^G aQ = 0, 
7?7f G af3 = 0. 

Setting A = det(X AB ) these equations can be put into the following form in the coordinates 

(x A ,x a ): 

_^ X l/2 D a {x -l/2 DaXAB) 

+ Xab( ~ \n + \- l ' 2 D a D a \ l l 2 + ^X-^DaX 1 / 2 - \D a \ CD D a \ CD ) = 0, 
Z 2 o 

(2.10) 

^D c X AB D d X AB -D c (\- l / 2 D d X 1 / 2 ) 

+ lcd {X- l l 2 D a D a X 1 ' 2 + l -D a X- l / 2 D a X 1 ' 2 - l -D a X AB D a X AB ) = 0. 

(2.11) 

Here D a is the derivative operator on M compatible with 7 a 5 and 1Z is the scalar curvature 
of D a . Indices for tensor fields on M (Latin lowercase) are lowered and raised with 7 a 5 
and its inverse r ) ab . It will be useful later to note that eq. (2.11) implies 

D a D a (X^) = 0. (2.12) 
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All of our considerations will be formulated on the space of orbits M. The Einstein 
equations (2.10), (2.11) are viewed as a generally covariant system of partial differential 
equations for "matter fields" \ab an d Lorentz metric 7 a ^ on the two-dimensional manifold 
M. 

2a. Cylindrical symmetry. 

Our model for open universes is defined by taking M. = R^ xS 1 and Q = S 1 x R. Note 
that M. is diffeomorphic to the manifold obtained by removing a timelike 2-plane (swept 
out by, say, the z-axis) from R 4 . The Killing coordinates are denoted by x A = ((f>,z), 
where <fi G (0, 27r) and z G (—00,00). The space of orbits is M = R x R+, where R + 
is the manifold of positive-definite real numbers. We define R := A2 and demand that 
V a R is everywhere spacelike. We will call spacetimes with these properties cylindrically 
symmetric. In what follows it will be useful to employ coordinates on M that are adapted to 
a foliation by spacelike curves R + M. These coordinates will be denoted by x a = (£, r) 
where t G (—00,00) and r G (0, 00). The radial variable r will be restricted by the 
requirement that ^ > 0. R itself can serve as a radial coordinate on R + . Because of 
(2.12), R will be a harmonic coordinate when the Einstein equations hold. The harmonic 
conjugate T to R is defined by 

D a T = e a b D b R, (2.13) 

where e ab) is the volume form on M defined by 7^. The coordinates t = T and r = R form a 
conformal coordinate chart when the Einstein equations hold, that is, in these coordinates 
the metric on M is given by 

ds 2 = Q 2 (-dT 2 + dR 2 ), (2.14) 
where O = 0(T, R) is a positive-definite function. 

We parametrize A^_g in terms of scalar fields (R 7 ip,ip) via 

X AB dx A dx B = R 2 e-^d(j) 2 + e^(dz + ^d<P) 2 . (2.15) 

In preparation for the Hamiltonian formulation of the Einstein equations (2.10) and (2.11), 
we foliate M with spacelike curves t = const and parametrize the metric on M via 

lab dx a dx h = [-(N^) 2 + e^(N r ) 2 ]dt 2 + 2e^N r dtdr + e^dr 2 . (2.16) 

As tensor fields on the curves t = const, N 1 - is a scalar field called the "lapse function", 
N r is a vector field called the "shift vector" (equivalent to a density of weight minus one 
on R + ), and 7 is the natural logarithm of a covariant rank-2 tensor field (equivalent to 
the logarithm of a density of weight-two). The functions (R,ip,ip) are scalar fields on 
t = const. All the fields (N^, N r , 7, R, -0, ip) are functions of the coordinates (t,r) only; 
see the Appendix for the boundary and fall-off conditions we use for these fields. In terms 
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of this parametrization, the conformal factor O arising in the coordinates x a = (T, R) (see 
(2.14)) is given by 

n = e^-^' 2 . (2.17) 

It is worth noting some special cases of the parametrization given above. If the Killing 
vectors are each hypersurface orthogonal, it follows that ip = 0, i.e., the matrix \ab is 
diagonal, indicating orthogonality of the Killing vectors. This special case of cylindrical 
symmetry in which each is hypersurface orthogonal will be called whole cylindrical 
symmetry. The solutions to the field equations possessing whole cylindrical symmetry are 
the well-known Einstein-Rosen waves. Kuchaf assumes whole cylindrical symmetry in his 
analysis [7] and, with ip = 0, our parametrization of the spacetime metric is identical to 
his. When 

N 1 - = 1, N r = 0, 7 = 0, R = r, V = 0, i> = (2.18) 

the spacetime is flat and the metric is expressed in cylindrical coordinates. 

Let an overdot denote differentiation with respect to t. The field equations (2.10), 
(2.11) in the parametrization (2.15), (2.16) can be obtained from the following Hamiltonian 
form of the action: 

SIN 1 , N T , 7, R, ?/>, 7T 7 , TV R , 7T^, 7T~] = jf dt dr (tT 7 7 + TV R R + TV^Ip + TT~^j 



*2 

dtn, 

tl 



(2.19) 



where the Hamiltonian is 

"OO 



H 



/■OO , , 

2Ar ± (oo)(l-e^°°)/ 2 ) + y dr [N^TL^ + N r H r ) , (2.20) 



and we have denoted the limits of AT-'-(r) and j(r) as r — > oo by AT-'-(oo) and 7(00). The 
super-Hamiltonian and super-momentum are given by 



H ± : = e 



. - P (^-7)/2 



- tvjTVji + 2R" - i?V + - (i?^ /2 + iT Vj) 



(2.21) 



7Y r : = -27r 7 + 7r 7 7 ; + 7TR.R' + 7ry,V + 7r ^V ,/ ' (2.22) 

where a prime indicates a derivative with respect to r. For each t the momentum variables 
(7r 7 , 7tr, 7r^, 7r~) are scalar densities of weight-one on R + . The Hamiltonian action func- 
tional (2.19) can be obtained by (i) expressing the ADM action in terms of the parametriza- 
tion (2.15) and (2.16), (ii) integrating out the Killing coordinates z and <j> (the range of z 
should be made finite), (iii) dividing the result by the range of z and the range of (p. 
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Extremising the action (2.19) with respect to variations of 7r 7 , 7Tr, tt^, and 7r~ produces 

the definitions of the momenta in terms of the "velocities" 7, R, ip, and ip. Extremising the 
action (2.19) with respect to variations of iV"- 1 -, N r , 7, R, ip, and ip, where the asymptotic 
values of N 1 - and N r are held fixed, yields six equations which, given the definitions of 
the momenta, are equivalent to the six field equations (2.10), (2.11) in the parametrization 
(2.15), (2.16). In particular, the field equations arising from varying the lapse and shift 
are the constraints 

H± « and H r ~ 0, (2.23) 

which correspond to the normal-normal and normal-tangential projections of the field 
equations (2.11) relative to the curve t = const. The remaining four equations are evolution 
equations for the spatial metric variable 7 and the "matter fields" R, ip, and ip. 

Let us make two remarks on the Hamiltonian variational principle we have described. 

(i) It is important to note that the canonical variables are restricted by the assumption 
that V a R is spacelike on M, which implies that D a R is spacelike on M. In Hamiltonian 
form, this restriction is 

R' > |tt 7 |. (2.24) 

(ii) The boundary term at infinity that appears in the Hamiltonian is needed to render the 
action and Hamiltonian functionally differentiable with the boundary conditions given in 
the Appendix. On solutions to the field equations the Hamiltonian is given by the boundary 
term, which we identify as the energy generating time evolution at infinity characterized 
by iV- L (oo). Note that our boundary conditions are such that iV- L (oo) is fixed, that is, not 
subject to variation. This allows us to add constant multiples of iV-'-(oo) to the surface 
term without spoiling the differentiability of the action or Hamiltonian. We have used this 
freedom to normalize the Hamiltonian so that it vanishes when spacetime is flat, which 
occurs when 7(00) = 0. The energy associated with time translations described by unit 
lapse function at infinity is 

£ = 2(l-e-^°°)/ 2 ). (2.25) 

It follows from the field equations that 7(00), known as the "C-energy", is conserved 
and non-negative (see, e.g., [19]). Hence E is bounded from below and is minimized on 
flat spacetime. Note also that E is bounded from above, which is due to the geometric 
interpretation of E as a deficit angle (divided by n) at infinity. See [20] for a detailed 
discussion of this interesting behavior of the energy of cylindrically symmetric spacetimes. 

To summarize, cylindrically symmetric spacetimes constitute a constrained Hamilto- 
nian system. The phase space V is the space of fields (7, R, ip, ip, 7r 7 , 7Tr, 7r^, 7r~) on R + 
with boundary conditions as specified in the Appendix, and with the restriction (2.24). 
The action functional (2.19) defines the symplectic 2-form O on T. O maps a pair of 



7 



tangent vectors 

X = (£7, SR, Sif}, Sip, Snj, 5n R , 5n^, 5n~) (2.26) 

and 

X = (?7, ?7T 7 , ?7Tfl, ?7T^, ?7T~) (2.27) 

to the real number 

tt(X, X) = J dr (Sir^Sj + Sn R 5R + Sir^Sip + 5tt~5$ - [5 < — ► 5]j . (2.28) 

From (2.28) it follows that (7, R, ip, ip) and (7r 7 , tt r , tt^,, 7r~) are, respectively, canonical 
coordinates and momenta for V and hence satisfy the canonical Poisson bracket relations, 
e.g., 

{ 1 (r),7T 7 (f)} = 8(r,f). (2.29) 

In terms of the Poisson bracket {•,■}, the time evolution of a function F: V — > R is given 
by 

F = {F,H}. (2.30) 

For any choice of lapse and shift, dynamical evolution takes place entirely on the constraint 
surface r T defined by the constraints (2.23). This follows from the fact that the 
constraints are "first class" . More precisely, the Poisson algebra of the super-Hamiltonian 
and super-momentum is the algebra of deformations of spacelike curves in M [21]. 

In two-dimensional spacetimes, such as we have here, it is convenient to work with a 
rescaled super-Hamiltonian which is a scalar density of weight-two. We define the weight- 
two super-Hamiltonian via 

H = eb-M 2 H ± . (2.31) 

Of course the super-Hamiltonian constraint H± ^ is equivalent to H « 0, and we can 
write 

N- L H L = NH, (2.32) 

where N is a scalar density of weight —1, which is equivalent to a vector in one dimension, 
obtained by 

N = e^-^/ 2 N L . (2.33) 



We can now vary the action 

TV. 1.. 7, 



S[N, N r , 7, R, Ip, 1p, 7T 7 , 7V R , 7ty, 7T; 



dt / dr (tt 7 7 + n R R + w^ip + tt-^ - NH - N r H^j - I dt JV (00)7(00), 

(2.34) 
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with respect to its arguments and obtain equations still equivalent to (2.10), (2.11) once the 
correspondence (2.33) is made. In this variational principle we hold fixed the asymptotic 
value of N. When using the weight-two super-Hamiltonian and associated lapse density 
the Hamiltonian is given by 

rOC 

H = AT(oo) 7 (oo) + dr (NH + N r H r ) , (2.35) 
JO 

where the surface term is again chosen to make S and H differentiable with our boundary 
conditions and to yield H = when 7(00) = 0. Evidently, the C-energy is associated with 
time translations defined by unit lapse density at infinity. 

The advantage of the weight-two super-Hamiltonian is that the Poisson algebra of 
constraints is now a Lie algebra [22]. In detail, define smeared constraints: 

roc poo 

H(N)= / drNH and H r (N r ) = / dr N r H r , (2.36) 

Jo Jo 

where N r is a given vector field and N is a given scalar density of weight —1. Direct 
computation then shows that 

{H(N),H(M)} = H r (J r ), 
{H(N), H r (M r )} = H(K), (2.37) 
{Hr(N r ),H r (M r )} = H r (L r ), 

where 

J r = NM' — MA/"', K = NM rl — M r N f , L r = N r M r ' - M r N r ' . (2.38) 



2b. Toroidal symmetry. 

Our model for closed universes is obtained by setting M. = R + x T^ and choosing 

the isometry group to be Q = T 2 with coordinates y, z G (0, 2n). The space of orbits is 

a cylinder, M = R + x S . We are considering one of the Gowdy models [15], and we 

shall call spacetimes with the above properties toroidally symmetric. In these spacetimes 
1 

we define r := A2. It can be shown that toroidally symmetric vacuum spacetimes are 
foliated by spacelike surfaces whose leaves are defined by r = const [15,23]. Hence the 
gradient of r is timelike, which we will assume in all that follows. According to observers 
"at rest" relative to the slices r = const, the toroidally symmetric spacetimes expand 
forever from a "big bang" at r = 0. Coordinates on M that are adapted to a foliation 
by spacelike circles will be denoted (t, x), where t G (0, 00) and x G (—00,00) with the 
identification x ~ x + 2n. We will demand that a time coordinate t on M satisfies ^ > 0. 
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For any solution of the Einstein equations, r is a harmonic time coordinate. The harmonic 
conjugate to r is denoted by X, which satisfies 



D * X = i r e l D a b \n ' (2 ' 39) 

where ds a is the oriented line element on a spacelike circle. Because D a r is timelike, the 
denominator in (2.39) never vanishes. By virtue of (2.12), the integral in the denominator 
is independent of the choice of spacelike circle, i.e., is a constant of motion. This integral 
is introduced so that on any circle t = const we have that 

X(x + 2n) -X(x) = 2n. (2.40) 

With the identification X ~ X + 2w, the coordinates t = r, x = X are adapted to a 
spacelike foliation of M; in these coordinates the metric on M is of the form 

ds 2 = Q 2 (-dT 2 + dX 2 ). (2.41) 



The Hamiltonian formulation of toroidally symmetric spacetimes closely parallels that 
obtained for cylindrically symmetric spacetimes. Essentially, the toroidally symmetric case 
differs by the change of notation R <-> r and the fact that the Cauchy surfaces in M are 
now compact circles parametrized by x instead of non-compact half-lines parametrized by 
r. We parametrize Xab m terms of scalar fields (r, tp, ip) via 

\ AB dx A dx B = T 2 e^dy 2 + e^(dz + ^dy) 2 . (2.42) 

The metric on M is parametrized relative to an arbitrary foliation by spacelike circles 
exactly as in (2.16): 

lah dx a dx h = [-(N^ 2 + e^(N x ) 2 ]dt 2 + 2e^N x dtdx + e^dx 2 . (2.43) 
In the coordinates x a = (r, X) the conformal factor in (2.41) is 

fl = e^-^l 2 . (2.44) 



When ip = our parametrization corresponds to a "polarized Gowdy model" [16], 
which can be considered the closed universe analog of the Einstein-Rosen waves. When 

N 1 - = 1, N x = 0, 7 = 0, r = t, V = 0, i> = (2.45) 

the spacetime is flat. This spacetime can be obtained from 4-dimensional Minkowski 
spacetime M 4 as follows. Let (£, x, y, z) be inertial coordinates on M 4 . Denote by I + the 
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chronological future of the origin in M 4 . Let I denote the manifold obtained from I + by 
the identifications 

x~x + 2n, tanh _1 (^) ~ tanh _1 (^) + 2tt, z ~ z + 2tt. (2.46) 

t t 

The Minkowski metric rj a p projects to a flat metric on I. The mapping 0: Ai — > I defined 

t = t coshy, 



(2.47) 



y = tsinhy, 

is a diffeomorphism which identifies the metric on Ai defined by (2.45), 

ds 2 = -dt 2 + dx 2 + t 2 dy 2 + dz 2 , (2.48) 

with the flat metric induced on /. 

In terms of a spacelike foliation of M with adapted coordinates (£, x), the Hamiltonian 
form of the action is 

^[A/^, iV r , 7 , T, ^, 7T 7 , 7T T , 7T^, 7T~] = dt dx (iT^ + 7T T T + 7T^ + TT-^J 

' J ° (2.49) 



ti 

*2 



dtn, 



where the Hamiltonian is 

"2tt 
/0 

The super-Hamiltonian and super-momentum are given by 



r 2ir 

H = J dx (N ± H ± + N X H X ) . (2.50) 



H ± : = e 



(V-7)/2 



-7r 7 7r T + 2r" - r' 7 ' + i (r<// 2 + r^n 2 ) + \ (re^V- + r~ l e 2 ^ 

« 0, (2.51) 
7^ : = -27r 7 + 7r 7 7 ; + tv t t' + ir^ip' + n-ijj' 

«0, (2.52) 

where a prime denotes differentiation with respect to x. As in the cylindrically symmetric 
case, the action functional (2.49) can be obtained by expressing the usual ADM action in 
the field parametrization (2.42) and (2.43), integrating out the coordinates y and z, and 
dividing by the ranges of these coordinates. Extremising (2.49) with respect to variations 
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of its arguments leads to the Einstein field equations (2.10) and (2.11) in our chosen 
parametrization. 

To summarize, points in the phase space V of toroidally symmetric spacetimes are 
defined by the smooth tensor fields on the circle (7, r, ip, 7r 7 , 7t t , 7r^, 7r~). These fields are 
restricted by the requirement that V a T is timelike on M., which means D a r is timelike on 
M; in Hamiltonian form this requirement is 

7r 7 <-|r / |. (2.53) 

The action (2.49) defines the symplectic structure on T. The symplectic 2-form acting on 
a pair of tangent vectors to V is given by 

tt(X, X) = J dx ^tt 7 ?7 + Sty t 5t + S^Sip + 8tt^ - [6 < — ► ?]) , (2.54) 

so that, with R <-> r, the canonical coordinates and momenta are as before. The weight- 
two super-Hamiltonian is defined via (2.31), and the Poisson-algebraic properties of the 
constraint functions are as in the cylindrically symmetric case. Dynamical evolution takes 
place on the constraint surface T defined by the constraints (2.51) and (2.52), and is 
generated by the Hamiltonian (2.50). 



3. Canonical Transformations. 

In this section we will exhibit canonical transformations from slight extensions of the 
gravitational phase spaces of §2 to phase spaces for parametrized field theories on a fixed 
background spacetime. The extensions are needed because the gravitational phase space is 
not quite adequate to define embeddings of hypersurfaces into Ricci-flat spacetimes. This 
difficulty arises for the full theory in both closed universes [14] and open universes [13], and 
can be considered a "global problem of time" [1] . While it is possible to phrase all of our 
results directly on the four-dimensional spacetime manifold M., it is far more convenient 
to express our results on the space of orbits M, and we will present our analysis on this 
effective two-dimensional spacetime manifold. 

Before specializing to our two Killing vector models, let us outline the basic strategy 
for the full theory [6,1]. Denote the usual gravitational phase space variables by (qij,p lJ ) 
and the super-Hamiltonian and super-momentum by (H±, Hi). We seek a canonical trans- 
formation 

( qiJ ,p lJ )^(X a ,U a ,q A ,p A ), (3.1) 

where a = 0, 1, 2, 3 and A = 1, 2 such that, on solutions to the equations of motion and 
constraints, X a : E — > M represents a spacelike embedding of a Cauchy surface E into the 
spacetime manifold M. The transformation must allow the constraints w Hi to 
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be resolved for the momenta conjugate to the embeddings, i.e., in the new variables the 
constraints are equivalent to 



H a :=U a + h a (X a ,q A ,p A )^0. 



(3.2) 



This formulation of the canonical theory has the following interpretation. The gravitational 
variables encoded in X a are used to identify instants of time and points of space at which 
the true degrees of freedom (<? a ,Pa) are being measured. The embeddings are "pure 
gauge", i.e., arbitrary; their conjugate momenta are completely determined in terms of 
the embeddings and true degrees of freedom by the constraints (3.2). The constraint 
functions H a , when integrated against functions N a , generate the dynamical evolution 
of the true degrees of freedom (q a ,Pa) as the embedding they are on is deformed, via 
8X a = N a , through the Ricci-flat spacetime for which (q A ,p A ) are Cauchy data. The 
densities of weight-one h a (X a , q A ,p A ) represent the energy-momentum current of (q A ,p A ) 
through the hypersurface embedded by X a . Gravitational dynamics on phase space is 
thus cast into the form of a "parametrized field theory" on M. 

3a. Open universes: cylindrical symmetry. 

As we saw in §2a, the metric variable R and its harmonic conjugate T define a conformal 
coordinate chart on cylindrically symmetric spacetimes (reduced to M). We can thus 
define a spacelike curve R + <— * M by giving its parametric description (T(r), R(r)). In the 
canonical formalism we can therefore view the phase space variable R(r) as one part of an 
embedding of R + into M. To complete the definition of the embedding we must express 
T as a function on the phase space, which will be again denoted T(r). Using the pullback 
of (2.13) to a spacelike curve t = const and the Hamilton equations we find 



Because of the restriction (2.24), the variables (T(r), R(r)) define a spacelike embedding 
on solutions of the equations of motion. Unfortunately, given a point in the gravitational 
phase space T, the embedding is not uniquely specified because the gravitational data do 
not fix the value of the integration constant T(oo), which represents the asymptotic location 
of the embedded curve. This difficulty is not a consequence of our use of the conformal 
coordinates (T, R) to specify an embedding, but is a general feature of general relativity 
of open universes. For example, in [13] the Hamiltonian formulation of general relativity 
in asymptotically flat universes is considered, and it is shown that the geometrodynamical 



T' = 



(3.3) 



This equation can be integrated to give 




(3.4) 
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data must be supplemented by a finite number of additional degrees of freedom in order to 
determine a spacelike hypersurface. To implement this idea in the context of cylindrically 
symmetric spacetimes, we will extend the phase space V of §2a. We then make a canonical 
change of variables on the extended phase space in which certain degrees of freedom are 
identified as spacelike embeddings. The resulting Hamiltonian structure of the model is 
then interpreted as a generally covariant form of a theory of harmonic maps by a sequence 
of canonical transformations on the appropriate harmonic map phase space. We begin by 
finding the appropriate extended phase space for cylindrically symmetric spacetimes. 

From the definition of the lapse density N and the time T, or equivalently, from the 
equations of motion for T, it follows that the rate of change of T(oo) with respect to t is the 
asymptotic value of the lapse density iV(oo). This motivates the following construction. 
We introduce a new degree of freedom, Too, which is the time displayed by a clock at 
infinity that measures time T. The lapse density at infinity is expressed as 



If we insist on keeping fixed the asymptotic value of the lapse density, then we also keep 
fixed foo, and nothing is changed except notation. However, we can treat Too as a new 
dynamical variable which is to be varied in the action. This is the usual logic of the 
parametrization process as applied to the "point at infinity" . If we parametrize at infinity, 
the Hamiltonian form of the action (2.34) can be written as 



By adding the new variable Too to the Hamiltonian action principle we get additional 
equations that, however, are equivalent to the original equations. In detail, by varying Too 
we obtain conservation of the C-energy, 7(00) = 0, which already followed from the other 
equations of motion and so does not alter the content of the field equations. The other new 
equation comes from varying 7. Prior to parametrizing at infinity, the variation of 7 led to 
one of the field equations (2.10), (2.11), and a potential boundary equation was eliminated 
by the boundary term in the Hamiltonian. After parametrization, the boundary equation 
survives and yields equation (3.5), which recovers the desired definition of Too- 

To summarize, we can enlarge the phase space Y of cylindrically symmetric spacetimes 
by adding a single variable Too- The extended phase space will be denoted r*. The extrema 
of the action functional (3.6) still define cylindrically symmetric vacuum spacetimes, but 
now in terms of the extended set of variables. By enlarging the phase space in this manner 
we are able to define the asymptotic location of spatial curves using dynamical variables. 
The action S, however, is not in Hamiltonian form because the surface term now enters as 



iV(oo) = foo- 



(3.5) 



S[N, N r 7 7, R, ?/>, 1p, 7T 7 , 7Tr, 7ty, 7T^, Too] 




(3.6) 
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a "kinetic term" and destroys the canonical nature of the phase space coordinates 7 and 
7TVy. We still must find canonical coordinates and momenta on r*. Indeed, we must show 
that r* is a symplectic manifold. We will take care of these issues, while at the same time 
providing the cylindrically symmetric version of the canonical transformation (3.1), in the 
following. 

Let us define the phase space T for a field theory on M as follows. A point in phase 
space is defined by the functions (T, R, tp, ip, IFj, H R , tt^, tt~) on R + , where (T,R,ift,ij)) 
are scalar functions and (IFf, H R , tt^, tt~) are scalar densities of weight-one. Boundary 
and fall-off conditions on these scalar fields and scalar densities are as indicated in the 
Appendix. We will need some restrictions on the functions T(r) and R(r) so that they can 
be interpreted as spacelike embeddings of R + into M. For reasons which will be clearer 
in a moment, we demand 

R'>\T'\, (3.7) 

and we include this inequality in the definition of T. We define a symplectic 2-form, 
S(X, X) on T by its action on a pair of vectors 

X = (ST, 5R, 64, 6$, 5U T , SU R , 8^, 8tt~) (3.8) 

and 

X = (ST, SR, SU T , SU R , Sn^, Sn~) (3.9) 

at a point of T. The symplectic form is defined by 

/*oo 

E(X, X) = J dr (SU T ST + SU R SR + Sir^Sip + Stt~S$ - [6 < — ► . (3.10) 

For the moment, the symplectic manifold (T, S) is to be viewed as logically independent 
of the gravitational phase space, that is, T is being thought of as simply a space of func- 
tions (T, R, ifj, if;, Ily, n^>, 7T.0, 7r~) upon which we have defined a symplectic structure. The 

symplectic structure is defined so that (T, R,ip,ip) and (Ily, U R , n^, tv~) are, respectively, 
canonical coordinates and momenta. 

Now consider the following map from T* to T. As suggested by our notation, we will 
identify the variables R, ip, ifj, 7r^, tv~ in V* and T. The remaining portion of the map is 
defined by 

T(r)=r 00 - f dfn^f), (3.11) 
J 00 

n T = -V + \\n(R /2 - irh] ' , (3.12) 



Hi? = ttr + 



7 

/R' -TV 

In 



R' -TT^ 1 ' 



R' + 71. 



(3.13) 
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Note that this transformation is consistent with and 11^ being scalar densities of 
weight-one. This map can be inverted; the relevant formulas are 



Too = T(oo), (3.14) 

rr 

7(r) = ln^' 2 - T /2 ) - / dfU T (f), (3.15) 

Jo 

tt 7 = —T f , (3.16) 

ta (*±rii'. (3.i7) 



i?' — T" 



The transformation (3.11)-(3.13) is a diffeomorphism which identifies and T. In par- 
ticular, the inequality (3.7) is precisely the restriction (2.24) on r*, and is needed for the 
transformation to be well-defined. We can use this diffeomorphism to express the action 
(2.34) as a functional S* of curves in T: 

S*[N, N r , T, Rrf, ^ U T , U R , tt v „ tt~] 

fh roo , . . (3.18) 

= J dt J dr \U T T + U R R + + tt-i/j - NH - N r H r ) , 

where 

H = U T R' + U R T ! + X - (r^' 2 + R~ 1 kI) + \ ^Re-^Tvl + IT V^' 2 ^) « 0,(3.19) 
H r = U T T' + T1 R R' + it^' + ix~i)' « 0. (3.20) 

Note that the surface term contribution to the action has dropped out. Indeed, from this 
action we see that the variables (T, R, ip, iji) and (Ily, II^>, 7r^, n~) are canonical coordinates 
and momenta for the phase space r*. We have (i) shown that T* is a symplectic manifold by 
exhibiting a diffeomorphism from V* to the symplectic manifold (T, S), and (ii) exhibited 
a canonical coordinate chart (T, R, ip, ip, Hy, IL R , tt^, tv~) on r*. 

The action S* has a nice mathematical interpretation in terms of a parametrized field 
theory formulation of harmonic maps on a flat spacetime, and we shall now spend a little 
time developing this interpretation. Recall that harmonic maps tp^: M — > L are fields on a 
spacetime (M, g) taking values in a Riemannian manifold (L, a) and which extremise the 
"energy integral" 

I[<P A ] = ~\! ^~99 a ^AB(^U% (3-21) 

where g a p is the metric on M, and a j± R is the metric on L. Note that we are using capital 
Latin indices to label the harmonic maps; in this discussion these indices should not be 
confused with those labeling the Killing vectors used to define the gravitational model. 
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Also, the spacetime manifold M and metric g a p for the harmonic map theory should not 
be confused with the gravitational M. and g a p. 

Introduce a spacelike foliation I:RxE^M characterized by lapse N 1 - and shift N\ 
The Hamiltonian form of the action (3.21) is 

I[<P A , *A] = [ Ua¥> A ~ N ± h ± - N'h,) , (3.22) 

where lowercase Latin indices denote tensors on E, and the energy and momentum densities 
are 

h± = \ ^cr AB n A n B + ^a AB ^fj , (3.23) 

hi = TTAtpj. (3.24) 

Here q^j is the induced metric on each hypersurface E of the foliation. Variation of 
I[(p^, t^a] with respect to cp A and -k a yields equations equivalent to those obtained by 
varying (3.21). At this point the spacetime metric and foliation, while arbitrary, are fixed. 
This is reflected by the fact that the lapse and shift are fixed fields on R x E, i.e., not 
subject to variation in the action principle. Because of this, the field theory — in either the 
Lagrangian or Hamiltonian formulation — is not "generally covariant" . General covariance 
can be introduced into the field theory by keeping the spacetime metric fixed and adding 
the foliation itself to the space of dependent variables to be varied in the action princi- 
ple. This is conveniently done in the Hamiltonian formulation, where the new dynamical 
variables are spacelike embeddings, which we shall denote by X a . We must still introduce 
momenta Ii a conjugate to the embeddings. To do this we need the unit normal n a to the 
hypersurface embedded by X a . The unit normal is defined by 

g aP (X)n a xP = and g a p{X)n a n^ = -1, (3.25) 

where g a p(X) is the metric on M restricted to the embedding X a . The unit normal is 
a fixed local function of the embedding and its first spatial derivatives. The foliation is 
a one-parameter family of spacelike embeddings, which we shall denote by X a (t). It is 
straightforward to show that 

X a (t) = N L n a + N l X a t . (3.26) 
The action (3.22) can thus be written as 

*A] = [ Ua<P A - X a h a ) , (3.27) 

where 

h a = ~n a h ± +X l a h h (3.28) 
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and we have introduced fields X a , which are fixed local functions of the embeddings and 
their first spatial derivatives defined by 

X* a X°j = 8) and n a X l a = 0. (3.29) 

From this form of the action we see that the momenta conjugate to the embeddings are 
given by 

n a = -h a ; (3.30) 

these definitions represent constraints 

H a := U a + h a ^ 0. (3.31) 

We can take the constraints (3.31) into account with Lagrange multipliers N a and obtain 
the final form for the Hamiltonian action describing the parametrized field theory: 

I[N a , X a , U a , </, n A ] = f (tt a <p A + U a X a - N a H a ) . (3.32) 

The extrema of this action, obtained by varying it with respect to its arguments, are 
defined by a system of equations equivalent to those obtained by extremising (3.27). Note 
that the Hamiltonian can be expressed as 

H := J N a H a = J (n ± H ± + N l H^j , (3.33) 

where 

H ± = n a Tl a + h ± , (3.34) 

Hi = X^Ua + hi, (3.35) 

and we can equally well vary N 1 - and N l instead of N a in (3.32). 

The resulting formalism, in which the embeddings and their conjugate momenta are 
adjoined to the phase space (at the expense of the constraints (3.31)) is the "parametrized 
formalism" for the harmonic map field theory. The central feature of the parametrized 
formalism is that it provides a generally covariant formulation of any field theory. This 
is what makes possible the identification of the gravitational models, which are generally 
covariant field theories on M, with the theory of harmonic maps on a fixed background 
spacetime. 

Now we are ready to make contact with the Hamiltonian formulation of cylindrically 
symmetric spacetimes. This is accomplished in 3 steps. 

1. Fix the spacetime and target space. 
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We fix the spacetime to be M = R 2 x S 1 equipped with a flat metric, which is defined 
in polar coordinates by the line element: 

g a(3 dx a dx (3 = -dT 2 + dR 2 + R 2 d<S> 2 . (3.36) 

Here, T G (—00, 00), R G (0, 00) and $ G (0, 2ir). The harmonic maps are defined to take 
values in L = R 2 with the metric <jab chosen to be of constant negative curvature. In 
coordinates ip G (—00,00) and 4> G (—00,00) the metric on L is determined by the line 
element 

a AB dip A dip B = di\) 2 + e~ 2i, d(f) 2 . (3.37) 
The scalar curvature of this metric is —2. 

2. Impose azimuthal symmetry. 

Of course ^ is a Killing vector field of the metric g a p. We now demand that the fields 
ip A be invariant along the flow generated by this Killing vector field, i.e., 

Ce_(p A = 0. (3.38) 
In the coordinates (T, R, $) on M this means we assume the fields are independent of 

ip = ij(T,R) and <f> = <f>(T,R). (3.39) 

In the parametrized formalism for the field theory we also assume that the foliation is 
compatible with the azimuthal symmetry in the sense that we only consider spacelike 
surfaces to which 4^ is everywhere tangent. The embeddings can be registered in the 
coordinates (T, R, $); that is, we have 

X a (r, $) = (T(r),R(r),$), (3.40) 

where r G (0, 00). Thus, to specify an embedding we must specify two functions of one 
variable, T(r) and R(r). On the symmetry-compatible foliation the shift vector takes the 
form 

AT* = (AT,0). (3.41) 

We can now formulate the parametrized harmonic map field theory on the 2-dimensional 
space of orbits M = R x R+ of the Killing vector. By working in polar coordinates on M, 
this amounts to simply ignoring the $ coordinate. The Hamiltonian form of the action can 
be obtained by substituting the choices made in the first two steps into the action (3.32) 
and integrating out the angular coordinate. We obtain 

/[JV-L, N r , T, RJl T , U R , V, 0, tt^, 7r ] 

t 2 roo 



rt2 roo , , 
= / dt dr (n T T + U R R + yr^ + yr^ - N ± H ± - N r H r ) , 

(3.42) 
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where 



H 



_L 



VR' 2 - T' 2 



U T R' + U R T' + l - [R^' 2 + iT Vj) 
+ \ (Re-^<P' 2 + R^n 2 ) 



(3.43) 



H r = U T T' + U R R' + n^' + (3.44) 

Variation of this action with respect to its arguments yields field equations equivalent 
to those obtained from (3.32) in the special case of azimuthal symmetry. An equivalent 
variational principle is obtained by defining a lapse density 

N := - = t N ± , (3.45) 

VR' 2 - T' 2 K ' 

and weight-two super-Hamiltonian 

H = U T R' + U R T' + \ (V 2 + R^n 2 ) + \ (ite^V 2 + R'^tt 2 ) . (3.46) 

In terms of these quantities we have 

/[AT, AT, T, R,U T , U R , V, (/>, 7r ] 

ft* f 00 / . x (3.47) 

= / dt dr \U T T + U R R + tt v > + tt^ -NH- N r H r ) , 
J t\ J 

and variation of this action with respect to its argument yields field equations equivalent 
to those obtained from (3.42). Note the close similarity between the super-Hamiltonian 
(3.46) for the harmonic map theory and the super-Hamiltonian (3.19) for cylindrically 
symmetric spacetimes. 

3. Canonical transformation. 

In the last step, we perform a canonical transformation that interchanges the roles of 
4> and 7T0 and puts the super-Hamiltonian and super-momentum of the parametrized field 
theory of harmonic maps into the form (3.19) and (3.20). This transformation is given by 

^(r)= f dfTrJf) (3.48) 



7T- = 0'. (3.49) 



The inverse transformation is 



r 



0(r) = J dfn~(f) (3.50) 

= (3-51) 
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With the boundary conditions given in the Appendix, this is a canonical transformation 
provided we impose the boundary condition <fi(r = 0) = 0. It follows immediately from 
(3.50) and (3.51) that, in the new variables, the super-Hamiltonian and super-momentum 
for the parametrized harmonic map field theory are precisely (3.19) and (3.20). 

As indicated by the equivalence to a parametrized field theory, the constraints (3.19) 
and (3.20), 

H ^ ^ H r , (3.52) 

can be solved for the momenta conjugate to the embeddings. Indeed, it is straightforward 
to show that these constraints are equivalent to 



H a : = Tl a + ha~ 0, 
where II a = (Hy, IIr), and h a = (/iy, hji) is defined as in (3.28) with 

R' T> 

n a = (n T , n R ) = 



VR r2 -T> 2 ' VR' 2 - T' 2 J ' 
T' R' 



R ,2 -T 2 ' R ,2 -T 2 
Explicitly, the constraints (3.53) take the form 

1 



(3.53) 

(3.54) 
(3.55) 



R i2 _ T /2 
- T' (tt^' + tt^') 



R> (i (V 2 + R-\l) + \ (Re- 2 ^ 2 - + R-'e 2 H' 2 ) 
*0, 



(3.56) 



and 



R' 2 - T' 2 
R' (tt^' + ir^' 



T> (\ (Rip* 2 + R-^l) + \ (Re- 2 ^l + R-'e 2 ^ 2 ) 



(3.57) 



0. 



A standard argument (see, e.g., [24]) establishes that the constraints (3.56) and (3.57) 
have an Abelian Poisson bracket algebra: 



{H a (r),H b (f)} = 0. 



(3.58) 



To summarize, the Hamiltonian structure of cylindrically symmetric gravitational fields 
is mathematically identical to a parametrized field theory of azimuthally symmetric har- 
monic maps from a 3-dimensional flat spacetime to a 2-dimensional manifold equipped 
with a metric of constant negative curvature. It is important to note that all of these 
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results are fully gauge invariant in the sense that no coordinate conditions are needed to 
be imposed on the gravitational theory. In any of the forms that we have presented the 
Hamiltonian formulation, the field theory retains the full 2-dimensional diffeomorphism 
invariance compatible with the imposition of cylindrical symmetry. 

3b. Closed universes: Toroidal Symmetry. 

We now repeat the analysis of §3a under the assumption of toroidal symmetry. The 
procedure is very similar to that used in the cylindrical symmetry case; the key difference 
is the way in which the missing degree of freedom is introduced. 

Again, the strategy is to use the conformal coordinates r and X to define embeddings 
of a circle into M = R + x S 1 . The variable r already appears as a canonical coordinate 
on the phase space, but we must still express X as a function on phase space. This can 
be achieved starting from (2.39). Choosing the spacelike circle on which the integral is 
performed to be a t = const slice, we have that 

X' = — tt 7 , (3.59) 

7TQ 



where 



1 



2tt 



7TQ := — / dxnj(x). (3.60) 
2tt Jq 

Note that, because of (2.53), 7r 7 is negative-definite, so the denominator never vanishes in 
(3.59) and X(x) is monotonic. 

Before completing the phase space definition of r and X into a canonical transfor- 
mation we have to contend with the fact that (3.59) does not define X as a function on 
the gravitational phase space V. The reason is the same as in the cylindrically symmetric 
case: there is an integration constant left unspecified in (3.59), which is not fixed by the 
gravitational phase space data. In §3a the integration constant represented the asymp- 
totic location of a spacelike slice, here it represents the relation between the origin of the 
coordinate a; on S 1 and the origin of the conformal coordinate X. This information is 
coordinate-dependent and not included in the gravitational phase space. As before, we 
remedy this situation by adding a new degree of freedom q to the phase space. Unlike the 
cylindrically symmetric case, this new degree of freedom has no role to play in gravitational 
dynamics, i.e., it is pure gauge. We therefore introduce a momentum p conjugate to q and 
adjoin a new constraint, 

p^O, (3.61) 

to the Hamiltonian formulation. We denote by V* the phase space extended by the variables 
q and p. Our strategy is to formulate the gravitational system as a parametrized field theory 
on r* and then reduce the system by the constraint (3.61). When we reduce T* by the 
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constraint (3.61) we arrive at the original gravitational phase space V and the dynamics 
thereon. 

We are now ready to extract embedding variables from the gravitational system for- 
mulated on r*. Let /ji(x) be a prescribed measure on the circle, i.e., fx is a positive density 
of weight-one on the circle, and is normalized via 



r2n 

/ dx /j,(x) = 1. 
Jo 



We define a transformation 



(T, 7, V', fa q, K T , 7T 7 , 7T^, 7T^, p) < ► (T, X, (p, <j>, Q, n T , U X , 7T^, V) , 



by 



7TQ 



/•27T /-a; | 

X(x) = q+ / dx" /j,(x") / dec' — ^(a/), 
JO A" ^0 

II T = -7TQ 7T T + 



In 



7r 7 — r 



n X = + ^0 ( l' ~ m (7l"7 - T ' 2 ) 



and 



= v 73 ^)^, 

4> = , — 
l 



TV. 



I r2n 
^0 Jo 



Q=ir I dx \ (t - Hrf ~ T ' 2 )) *i ~ ( + 
1 1 



In 



7r 7 — r 



r' + 7r. 



7 



P = 7T . 

ave mad 

the definition of X guarantees that 



3.62) 
3.63) 

3.64) 
3.65) 
3.66) 
3.67) 

3.68) 
3.69) 

3.70) 
3.71) 



3.72) 
3.73) 



Here we have made some convenient rescalings. In particular, the scaling by ^ used in 



X(2n)-X(0) = 2tv. 
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(3.74) 



We mention that the requirement that D a r is timelike, given by (2.53), is equivalent to 

X' > \T'\. (3.75) 

The inequalities (2.53) and (3.75) guarantee that the above transformation is well-defined 
and that the slices embedded by (T(x),X(x)) are spacelike. Let us also note that the 
constraint (3.61) is, in the new variables, the constraint 

r2ir 

/ dxU x {x)^0. (3.76) 

Jo 

It is not too hard to show, e.g., by expanding in Fourier series, that the above transforma- 
tion is a bijection from T* to the phase space T for a field theory of (0, 0, <2, 71^, tv^ V) in 
the parametrized formalism. Our notation here is that T is the product of the cotangent 
bundle over the space of embeddings of a circle into M, and the phase space of the canon- 
ical variables (</>,</>, Q) and {tt^tt^V). By computing Poisson brackets, or by computing 
the symplectic form in the new variables, it is straightforward to verify that the trans- 
formation is canonical, i.e., identifies the respective symplectic structures, and that the 
variables (T, X, <f>, <j>, Q) and (Ily, 11^, n^, V) are canonical coordinates and momenta. 

Modulo the constraint (3.76), the weight-2 super-Hamiltonian and super-momentum 
take the following form when expressed in terms of the new canonical variables: 

H := U T X' + U X T' + \ (V 2 + Iirj) + \ (Te-^Kl + » 0, (3.77) 

and 

H x := U T T' + Ii x X' + n^' + tt^' » 0. (3.78) 

The constraints of the theory are (3.76)-(3.78); they are "first-class". In particular, the 
Poisson brackets of the constraint (3.76) with the super-Hamiltonian and super-momentum 
vanish because (3.76) generates constant shifts of X with respect to which H and H x are 
invariant. 

Let us make one further remark about the canonical transformation (3.63). First, 
we emphasize that the identification of T* with T depends on the non-gravitational data 
(q,p,f/,). However, upon passing to V using the constraint (3.76) this dependence neces- 
sarily disappears. To see this, recall that the reduction from T* to V is obtained by (i) 
restricting the phase space to the constraint surface defined by (3.76), and (ii) identifying 
points on the constraint surface which lie on an orbit of the canonical transformations 
generated by J s i Tlx, that is, embeddings X and X + const are identified. The dependence 
of on p and [i is eliminated upon passing to the constraint surface defined by (3.76). 
Restricting to the constraint surface of (3.76), we must identify X and X + const. For a 
fixed choice of ji this eliminates the dependence of the phase space on q. Now suppose we 
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use a different measure fi to define a new embedding variable X(x). Because X' = tv^/tvq 
we have that (X — X)' = 0, that is, X = X + const, which is precisely the transformation 
generated by J s i Hx, so on T the phase space is independent of the choice of measure [i. 
Put differently, changing \x is equivalent to holding // fixed and changing q, which has no 
effect on points of V. Thus the non-gravitational data ((/,£>, /(/) are eliminated upon passing 
from r* to the original phase space V. 

Once again we can interpret the resulting formalism in terms of a parametrized field 
theory of harmonic maps. However, it must be kept in mind that there is a "point particle" 
degree of freedom represented by Q and V along with an extra constraint (3.76). Because 
Q is cyclic in the Hamiltonian, we can reduce the phase space V* by the integral of motion 
V to a phase space T*. On T* we view V as a parameter, which will appear in the metric 
for the spacetime upon which the harmonic maps are defined. To express the gravitational 
theory formulated on T* as a parametrized harmonic map field theory we repeat the 3 
steps of §3a. 

1. Fix the spacetime and target space. 

The spacetime is now taken to be M = R + x T 2 with metric given by the line element 

ds 1 = I 2 {-dT 2 + dX 2 + T 2 dY 2 ^j , (3.79) 

where X and Y are coordinates on T 2 , T > 0, and / is a positive constant. Using an 
analogous construction to that found in §2b, this metric can be viewed as a flat metric on 
a compactification of the chronological future I + of the origin of 3-dimensional Minkowski 
spacetime M 3 . Let us spell this out in detail. Let i, x, y denote inertial coordinates on 
M 3 . On I + make the identification 

x~ x + 2nl, tanh _1 (^) ~ tanh _1 (^) + 2tt, (3.80) 

t t 

and define the resulting manifold by I. The Minkowski metric on M 3 projects to a flat 
metric on /. Define coordinates (T, X, Y) on R + x T 2 , where T e (0, oo), X e (— oo, oo), 
and Y E (-oo, oo), with X ~ X+2ir and Y ~ Y+2n. The diffeomorphism 0: R+xT 2 ^ I, 
defined by 

i = IT cosh Y, 

x = lX, (3.81) 
y = ZTsinhy, 

identifies the metric (3.79) on R + x T 2 with the Minkowski metric on /. 

As in the cylindrically symmetric case, the harmonic maps cp^ := (ct,/3) take values in 
L = R 2 , which is equipped with a metric of constant negative curvature given by the line 
element 

a AB d^ A dp B = da 2 + e- 2a d[5 2 . (3.82) 
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2. Impose azimuthal symmetry. 

The vector field |y is a Killing vector field for the metric on M. We now demand that 
all the fields of the parametrized harmonic map field theory are likewise invariant under 
the f -parameter family of isometries generated by Sjy . Thus we can formulate the theory 
on the space of orbits, M = R + x S 1 , of the Killing vector field. On the space of orbits of 
?jY j the Hamiltonian form of the action is given by 

I[N, N x , T, X,n T , n x , a, 0, n a , n p ] 

ft* f 2lT / ■ ■ ■ \ (3.83) 
= J dt J dx [tv t T + tv x X + iv a a + Txp(3 - NH - N X H X J , 

where 

H = n T X ! + n x T ! + 1 ((/T)" 1 ^ + ITa' 2 ) + \ ((/T)"^ 2 ^ + lTe' 2a f3' 2 ) , (3.84) 
and 

H x = % T T' + n x X' + n a a + irpP . (3.85) 

Note the similarity between the super-Hamiltonian (3.84) for the harmonic map theory 
and the super-Hamiltonian (3.77) for toroidally symmetric spacetimes. 

3. Canonical transformation. 

In this last step we make a canonical transformation which puts the super-Hamiltonian 
and super-momentum (3.84) and (3.85) of the parametrized field theory into the form (3.77) 
and (3.78) found in the gravitational theory. It is possible to adapt the transformation 
(3.48)-(3.49) used for the cylindrically symmetric case, but this leads to additional— 
and unnecessary — constraints on the harmonic maps. In the toroidal symmetry case the 
following canonical transformation accomplishes our goal [25]. 

a = -r^ + ln(ZT), (3.86) 
= (3.87) 
ir a = -l^cji + lX', (3.88) 

Trp = Htt^, (3.89) 

n T = n T + l^T- 1 n^-l(2T)- 1 X', (3.90) 

n x = H x + l^4>' -l{2T)- l T' . (3.91) 

Note that the embedding coordinates (T, X) retain their original meaning, and the em- 
bedding momenta, while redefined, are still scalar densities of weight-one on the circle. By 
using the transformation (3.86)-(3.91) it follows that the super-Hamiltonian and super- 
momentum (3.84) and (3.85) of the parametrized field theory become those found in the 
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gravitational theory formulated on V* ((3.77) and (3.78)), provided we make the identfi- 
cation 

/ = -V. (3.92) 

As in the cylindrically symmetric case, we now know that the constraints (3.77) and 
(3.78), 

H&Q&H X , (3.93) 

can be reexpressed as 

H a := n a + h a ~ 0, (3.94) 
where n a = (Ily, IT^), and h a = (hj<, h x ) is defined as in (3.28) with 

X' T' 



n a = (n T ,n x ) 



^JX' 2 -T' 2 ^X' 2 -T' 2 
T' X' 



X>2 _ j>/2 ' x' 2 - T' 2 
Explicitly, the constraints (3.94) take the form 
1 



(3.95) 
(3.96) 



X' 2 - T' 2 



X' Q (T<p> 2 + T- 1 ^) + \ (Te- 2 ^/^? + T- 1 e 2 ^/^ /2 )) 



+ tti<; 



o, 



(3.97) 



and 



n 



X 



X' 2 - T' 2 



- X' [ixxd) + it id! 



0. 



(3.98) 



As before, the constraints (3.97) and (3.98) have an Abelian Poisson bracket algebra: 

{H a (x),H b (x)} = 0. (3.99) 

To summarize, the Hamiltonian formulation on V* of toroidally symmetric gravitational 
fields is equivalent to that of a parametrized field theory of harmonic maps from a flat 3- 
dimensional spacetime to a 2-dimensional manifold of constant negative curvature. Of 
course, to recover the original gravitational field theory on V the degrees of freedom (<2, V) 
must be reinstated but, more importantly, the constraint (3.76) must be imposed on the 
harmonic map field theory. This constraint sets to zero the homogeneous mode of the 
embedding momentum il^, and, by virtue of (3.98), is equivalent to the constraint 

"2tt 

dxh x (x)^0. (3.100) 
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This constraint can be viewed as setting to zero the total momentum of the fields <p and <f> 
(or a and f3). Because the constraint (3.100) is first-class, we are obliged to identify X and 
X + const. The equivalence class of embeddings thus obtained is independent of the value 
of the non-gravitational variables (q,p, (i). On the other hand, the extra constraint means 
that, strictly speaking, the Hamiltonian formulation of toroidally symmetric spacetimes is 
not quite identical to a parametrized field theory of harmonic maps. However, as we shall 
see, the correspondence is sufficiently close to allow a formal construction of the quantum 
theory. 

4. Implications for quantum theory. 

Issues which arise in canonical quantum gravity, the problems of time in particular [1] , 
can be investigated in the models we have been studying. In the last section we exhibited 
canonical transformations that identify certain phase spaces for these models with phase 
spaces for harmonic maps from flat 3-dimensional spacetimes into a target space of constant 
negative curvature. This mathematical identification is fully gauge invariant, i.e., preserves 
the 2-dimensional diffeomorphism symmetry exhibited by the two Killing vector models. 
It should be emphasized, however, that the extraction of the embedding variables was 
not without its complications. We encountered a "global problem of time" which was 
resolved by adding a finite number of non-gravitational degrees of freedom to the usual 
ADM phase space. Given the extended phase spaces, each event in the effective two- 
dimensional spacetime is uniquely labeled by the values of the canonical variables X a on 
a spacelike slice at the point where the spacelike slice passes through that event. The 
identification of spacetime points provided by X a is independent of the choice of slice 
and the "spacetime problem" [1] is avoided. Given these results, the obvious strategy 
for quantization is to view the dynamics of the models as a generally covariant (i.e., 
parametrized) formulation of non-linear fields on a fixed background. The quantum theory 
of fields on a fixed spacetime is quite well-studied and, at least superficially, presents no 
overwhelming conceptual difficulties — although the quantum theory of interacting fields 
is always technically challenging. A simplifying feature of these models is that, while 
they can be viewed as field theories on a 3-dimensional flat spacetime, the Killing vector 
structure is such that the final result is in each case a 2-dimensional field theory, and in two 
dimensions quantum field theory is typically more manageable than in higher dimensions. 
Our purpose in this section is to discuss certain broad features of the quantization of these 
models based upon the classical structures elucidated in §3. We hope to return to a more 
detailed examination of the resulting quantum theories in future work. 

In each of the models we have studied the resulting Hamiltonian structure involves 
embedding variables X a and their conjugate momenta n a , along with dynamical variables 
(</ a ,Pa) which are (ip, ip, 7r^, 7r^) in the cylindrical symmetry case, or (Q, 0, 0, V, tt^, irj) 
in the toroidal symmetry case. The constraints include diffeomorphism constraints (3.53), 
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(3.94) and, in the case of toroidal symmetry, an additional constraint (3.76). Ignoring 
its origins, the parametrized field theory defined by these constraints is, at least formally, 
relatively straightforward to "quantize". Note however that in each of the models the 
effective 2-dimensional field theory has a complicating feature not usually found in more 
familiar 2-dimensional field theories. Namely, the super-Hamiltonian is an explicit function 
of R(r) in the cylindrical symmetry case, and T(x) in the toroidal symmetry case. This 
will modify the quantum theories relative to what we might find, e.g., for a conformal 
field theory. In any case, we promote the canonical variables (q a ,Pa) to operators on 
states >. One way to do this would be, using the Heisenberg picture, to quantize the 
variables (q A ,p A ) on the priveleged foliation T = t and R = r or X = x. The states 
are embedding-independent in the Heisenberg picture; the operators are evolved by the 
constraint operators [26] . Of course it may be necessary to use perturbation theory to define 
the operators and states. Formally, observables are self-adjoint operators representing 
classical functionals of (q A ,PA)- Note that in the toroidal symmetry model, observables 
must commute with the constraint (3.76). 

In the Schrodinger picture, dynamical evolution corresponds to considering states 
|\&,X a > which are parametrized by the embeddings. The states are evolved from one 
embedding to the next by the energy- momentum current h a . More precisely, the states 
|W,X a > are defined as solutions to the functional Schrodinger equation 

^,X a >= h a \^,X a > . (4.1) 
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This equation can be considered an implementation of the demand that "physical states 
are annhilated by constraints": H a \^>,X a >= 0. In the toroidal symmetry case we must 
also impose the quantum version of (3.76): 

' dx lYM l*^ a >=0- (4-2) 
bX{x)J 

This requirement is equivalent to 

\V,T,X+ const >= \V,T,X >, (4.3) 
but, in light of (4.1) is also equivalent to 

(^j dx h x (x^j \y,X a > = 0. (4.4) 

Given an initial state \^q,Xq > on an m itial embedding Xq, the functional Schrodinger 
equation (4.1) is solved (subject to the subsidiary condition (4.2)) and the solution is 
matched to the initial data. Note that (4.2) can be implemented by imposing (4.3) on 
the initial state \^q,Xq >, and then solving (4.1). If (4.2) is satisfied initially, it will be 
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satisfied on any embedding provided the state is evolved according to (4.1). The resulting 
embedding-dependent vector X a > is to be interpreted as the state of the system on the 
embedding X a (x) and one can go on to predict outcomes of measurements of observables 
constructed from the quantum fields (and X a (x) if desired) on various hypersurfaces in 
the usual way. 

For this quantization strategy to be feasible, one must check that the operator represen- 
tatives of the currents h a are well defined operator-valued distributions. Furthermore, one 
must guarantee that the infinite number of differential equations contained in (4.1), (4.2) 
are mutually consistent. The integrability conditions for this are the quantum commutator 
analogs of the Abelian algebra (3.58) and (3.99), which guarantee that quantum dynamics 
does not depend on the choice of observer. More precisely, the integrability conditions for 
(4.1) imply that the state Xf > obtained by evolving an initial state \^>q, -^o ^ along 
a specific foliation connecting the slices defined by Xq and X^ does not depend on the 
choice of the foliation. For these integrability conditions to be satisfied, it is essential that 
the commutators of the components of the energy-momentum current are, up to a factor 
of z, the same as their classical Poisson brackets. It is very likely that this will not happen, 
i.e., anomalies (or "Schwinger terms") terms will arise. This is the "functional evolution 
problem" [1]. It would be very interesting to compute these Schwinger terms. If they 
depend only on the embedding, and are finite, then one can use techniques developed by 
Kuchaf for free fields in two dimensions, where the anomalous terms are finite and depend 
only on the embeddings [26]. If the anomalous terms are operator- valued, or infinite, it 
is not at all clear how to proceed. Note that if we turn off one of the polarizations of 
the gravitational field by setting ip = 0, then the resulting field theory is linear. In this 
case it would seem that all operators can be defined by normal ordering, and that the 
Schwinger terms are finite and depend only on the embeddings, but to our knowledge no 
one has computed the Schwinger terms for this interesting special case. In many ways, the 
functional evolution problem is the most important issue to address in studying the quan- 
tization of the two Killing vector models, for it is here that one verifies the compatibility 
of the method of quantization with the principle of general covariance. 

Given a successful resolution of the functional evolution problem, we should relate 
the elements of the quantum parametrized theory with geometrical elements of spacetime. 
In particular, what does it mean to "know the state on the embedding X a (x) nf ? Let 
us treat each of the models in turn. To fix the state |W,X a > of a cylindrically sym- 
metric spacetime, we imagine a family of observers and the associated spacelike foliation 
of spacetime. On a given slice of the foliation the observers measure the gravitational 
variables (R, 7r 7 ) and the "laboratory variable" Too; this fixes the embedding of the slice, 
X a (r) = (T(r), -R(r)), via (3.11). The observers also measure measure a complete set of 
commuting observables built from the operator representatives of (1^,1(^,7x^,7:^). All to- 
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gether, these measurements fix the quantum state on a spacelike hypersurface in terms 
of measurement of geometrical quantities and the reading of a clock at infinity. In the 
toroidal symmetry case the family of observers measure t/ttq and ti^/tyq; this determines 
the embedding X a (x) = (T(x), X(x)) via (3.64) and (3.65) up to specification of the non- 
gravitational variables q and //. The observers also measure a complete set of commuting 
observables built from (Q, cj), 0, V, n^,, tt^), which are defined in terms of the spacetime ge- 
ometry in (3.68)-(3.73). Different choices of q and fx simply redefine X(x) by addition of a 
constant. Because the states of a toroidally symmetric spacetime must be invariant under 
this transformation (see (4.3)), the quantum state \^,X a > is unambiguously determined 
by measurements of gravitational data only. 

There is one remaining problem that should be discussed when quantizing general 
relativity as a parametrized field theory. Kuchaf calls this problem the "multiple choice 
problem" [1]. The quantization we have outlined apparently depends quite heavily on the 
way in which the canonical variables X a are used to identify points of space and instants 
of time. To be sure, the embeddings we have constructed are, in some sense, geometrically 
natural. But it is conceivable that other, equally valid, embedding variables could be con- 
structed. In this case it is not at all clear that when using some other embedding variables 
the resulting quantum theory will coincide with the one we have outlined here. If the quan- 
tum theories based on different choices of embedding variables do not coincide, i.e., are 
not physically equivalent, then we have an embarassment of riches: which quantum theory 
describes the real world? This issue can be examined in the context of the parametrized 
formalism for a relativistic particle moving in a curved spacetime, and it is found that the 
multiple choice problem can be quite severe [1]. In many ways it is the multiple choice 
problem which most deeply reflects the conflict between general relativity and quantum 
mechanics. In order to examine this problem in the models we have considered, new sets of 
embedding variables are required. For example, one relatively simple way to obtain other 
embedding variables is to take the embeddings X a we have constructed here and perform 
the point transformation 

X a :=F a (X). (4.5) 

Assuming F a is smooth and admits a smooth inverse, this transformation can be easily 
completed to a canonical transformation. These new embeddings can be interpreted as 
follows. Fix an embedding X a of a hypersurface E in spacetime. Displace E by letting the 
diffeomorphism F a act on it pointwise. The resulting hypersurface is embedded by X a . Is 
it possible to construct a quantum theory of the parametrized field theories representing the 
classical two Killing vector models so that transformations such as (4.5) lead to physically 
equivalent theories? Of course we cannot answer this question here since we have only 
outlined the most basic features to be expected in the putative quantum theory. Let us 
however suggest that the multiple choice problem is closely allied with a familiar issue that 
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arises in quantization of a gauge theory using gauge fixing conditions. There, the question 
is whether or not the predictions of quantum field theories based on different gauge fixing 
conditions agree. The connection of this issue with the multiple choice problem stems from 
the fact that, at least in the models considered here, the classical gravitational theory can 
be obtained by (i) choosing the gauge T = t and R = r or X = x, and (ii) "parametrizing" 
the resulting field theory to regain general covariance. Likewise, we quantize the fields 
(q A iPa) on the foliation T = t and R = r or X = x, and then reinstate general covariance 
by constructing states parametrized by the embeddings X a which satisfy (4.1). Whether 
general covariance is in fact realized in the quantum theory is determined by the outcome 
of the problem of functional evolution. Let us assume that this problem can be solved. 
The remaining question is how the quantum theory depends on the initial choice of gauge. 
In gauge theories, this issue is fruitfully analyzed using BRST methods [27] which, at least 
formally, guarantee the independence of the predictions of the quantum theory from the 
underlying choice of gauge. It remains to be seen whether the multiple choice problem in 
generally covariant theories can be resolved in a similar way. 
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Appendix: Boundary conditions for cylindrically symmetric spacetimes. 

In this Appendix we summarize the boundary conditions used in the analysis of cylin- 
drically symmetric spacetimes. These boundary conditions guarantee existence and differ- 
entiability of the various action and Hamiltonian functionals that are used in the paper, 
as well as the canonical nature of the various transformations considered. The boundary 
conditions are also such that they are preserved under the dynamical evolution generated 
by the Hamiltonian functionals. It is conceivable that weaker boundary conditions can be 
used, but those used here allow for a reasonably large class of solutions to the field equa- 
tions. In particular, the boundary conditions include the Einstein-Rosen wave solutions 
which arise when one assumes whole cylindrical symmetry. 

It is important to note that the asymptotic values (r — > oo) of the lapse and lapse 
density, while arbitrary, are held fixed when varying the Hamiltonian (or action) defined 
on the ADM phase space V. The asymptotic values of the lapse and lapse density are 
only allowed to vary when using the extended phase space T* for cylindrically symmetric 
spacetimes. Note also that many other variables have a non-vanishing value as r — > or 
r — > oo; these limiting values are not held fixed in any of the variational principles. In 
particular, the values of the lapse and lapse density on the axis r = are not held fixed. 

The phase space variables are functions of the radial coordinate r; we must give the 
boundary conditions as r — > and r — > oo. As r — > we assume the following behavior of 
the canonical coordinates and momenta on T: 



R 



7 



r + C(r 3 ), 
^(0) + O(r 2 ), 
O(r), 



(Al) 
(A.2) 
(A3) 
(A.4) 
(A.5) 
(A.6) 
(A.7) 
(A.8) 






The lapse and shift are assumed to have the following behavior as r — > 0: 



= N ± (0) + O(r 2 ), 
N r = 0(r 2 ). 



(A.9) 
(A.10) 



The lapse density inherits its behavior as r — > from the lapse function: 



N = N(0) + O(r 2 ). 



(All) 
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Equations (A. 6) and (3.4) imply that as r — > 0, 

T = T(0) + C(r 3 ). (A.12) 

Similarly we assume 

n r = 0{r) (A.13) 

and 

n# = 0(r). (A.14) 

Our boundary conditions at r = imply that there are no singularities on the axis of 
symmetry. 

As r — > oo we assume the following behavior of the canonical coordinates and momenta 
on T: 



R = r + (D(r- e ), 


(A.15) 


7 = 7(00) + 0(r~ e ), 


(A.16) 




(A.17) 


$ = 0(r- e ), 


(A.18) 


n R = 0(r- 1 ), 


(A.19) 


7r 7 = 0(r-( 1+£ )), 


(A.20) 


n^ = 0(r- 1 ), 


(A.21) 




(A.22) 



where e > 0. The lapse and shift are assumed to have the following behavior as r — > 00: 

N L = JV-^oo) + 0(r- e ), (A.23) 

jV r = 0(r- £ ). (A.24) 

The lapse density inherits its behavior as r — > 00 from the lapse function: 

iV = A^(oo) + C(r" e ). (A.25) 

Equations (A.20) and (3.4) imply that as r — > 00, 

T = T(oo) +0(r- e ). (A.26) 

Similarly, 

n r = 0(r-( 1+e )) (A.27) 

and 

= 0{r- 1 ). (A.28) 
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